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The parton content of the electron is analyzed within perturbative QCD. 
It is shown that electron acquires an anomalous component from QCD, anal- 
ogously to photon. The evolution equations for the 'exclusive' and 'inclu- 
sive' electron structure function are constructed and solved numerically in 
the asymptotic region. 

I. INTRODUCTION 

The photon structure function describes the distribution of QCD partons inside a photon. 
It is known for long [|l| to have 'anomalous' component, which is calculable within pertur- 
bative QCD and dominates at asymptotically large momentum scales. This asymptotic 
solution, as opposed to those for hadrons, is independent of input data measured at lower 
momentum scales. At finite scales the photon structure gets modified by both perturbative 
and non-perturbative QCD contributions. 

The QCD structure of the photon is revealed in interactions with a highly virtual 'probe'. 
To fix attention let us think of a virtual gluon G* with momentum q, probing the photon 
which gets resolved into QCD partons. Their density f^{x,Q^) {k = q,q,G), depends on 
fractional momentum x of the parton with respect to photon and on the gluon virtuality 

= which must be large as compared to the QCD scale Aq^j-,. 

The photon structure is measured in experiments where the electron serves as a target. 
The process is depicted in Fig.|l|a, where also the notation is given. The black blob denotes 
'resolved' photon and sums up all coUinear QCD contributions. The full cross-section gets 
also a contribution from the hard ("direct") G*'j scattering (see e.g. 0), but we will not 
discuss it in this paper. Actually the photons emitted by the electron are virtual and, from 
the point of view of a physical process, G* measures the structure (parton content) of the 
electron, as depicted in Fig.p^. 





e(p') ^^^^ e(p') 

FIG. 1. Deep inelastic scattering on a photon (a) and electron (b) target 
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The aim of this paper is to study the QCD predictions for the electron structure function 
at large momentum scales. In particular we will discuss the dependence on the maximal 
virtuality of intermediate photons emitted by the electron, i.e. on the maximal momentum 
transfer between final and initial electron. In the next section we present the problem in 
general and discuss the relation to experimentally measured quantities. In Sec. 3 we derive 
the QCD solution to the non-singlet electron structure function in the moments space. A 
complete set of evolution equations is constructed and solved in Sec. 4. In Sec. 5 we present 
the summary and outlook. 



II. GENERAL FRAMEWORK 

The space-like virtuality of the photon exchanged in the diagram in Fig. |I]a 

(p-pT^-p^ (1) 

can be fixed by measuring the momentum p' of the outgoing electron. Otherwise it lies 
within kinematic limits 

PLM - < < Q^ ' + y -'y , (2) 

i. y z 

where me is the electron mass, y = qp-y/qp and z = Q^/2pq (these are approximate expres- 
sions valid at — the exact formulae can be found in 0). Usually the upper limit 
on P^ is set by experimental conditions (e.g. by anti-tagging) 

< P'. (3) 

The density of partons {k = q, q, G) with momentum pk = zp seen by our probe in the 
electron reads 

p2 

ft{z,Q\P^) = J dxdy6{z-xy) J dP^ f;{y, P^)f,{x,Q\ P^) (4a) 

P'^- (y) 

mill ■' 

= /v Jdp'f;iy^p')ml,Q^p',), ^^^^ 



V 

mm ' 



where 



27r P2 



y Pf 



(5) 



and f^{x,Q^, P^) describes the (j*7 interaction. In Eq.(^ only transverse photons are 
taken into account which is correct within the leading order of perturbative QCD. All QCD 
contributions are contained in f^{x,Q^, P^), which is discussed in the literature as the 
structure function of virtual photon @J^J^. In the following we will assume that P^ ^ 
which allows us to neglect the second term in the square brackets of Eq. (Kl) . 
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For Q ^ P f2 (x, Q , P^) can be approximated by the structure function of real photon 

)2 

7 
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[P^ = 0) and upon integration over we arrive at the Weizsacker-Wilhams formula: 

dy let ^2 r.^ ^ P 



ft{z, Q\ p') ^ f -^f^{y)m^, g^ 0) log (6) 



where 



f-^y^ = ^ — -y — ■ 

Formula has probabilistic interpretation in terms of the density of photons emitted by 
the electron and the density of QCD partons within the photon. As discussed in this 
partonic picture breaks down at very high energies when Z and W bosons contribute. 

In general, the experimentally measured electron structure function is always integrated 
over a range of photon virtualities and summed over the contributions from all weak in- 
termediate bosons. This structure function describes the QCD content of a real (on-shell) 
electron and allows for probabilistic interpretation of the cross sections. Even when we ne- 
glect the contributions from Z and W bosons the integration over P^ disables the partonic 
interpretation [|^. As compared to the standard QCD structure functions the electron one 
has extra dependence on maximal photon virtuality P^, which means that we do not inte- 
grate over all final electron states. In this sense we say that this electron structure function 
is 'exclusive'. 

Eq.(|^) is an approximation to the electron structure function for ^ P^. The approach 
to this limit within QCD is discussed in the next section. 

In the following we will assume that the Z and W bosons do not contribute but we 
will allow for arbitrary P^ < Q^. For both and P^ much greater than we have (c./. 
Eq.(i) 

p2 

m\p')= I ^f^^fAQ^p^), (8) 

.2 ^7 



mill 



with explicit z dependence suppressed and ® denoting convolution 



1 1 

(/ ® 9)iz) = dx dy6{z - xy) f{x)g{y). (9) 







We know from experiment that a nearly real photon (P^ ^ ^qcd) has a hadronic compo- 
nent which is often described phenomenologically in terms of the Vector Meson Dominance 
model (VDM) (see e.g. and references therein). This non-perturbative hadronic com- 
ponent becomes less important at higher Q^. As will be shown in the next section, any 
perturbative QCD predictions for P^ dependence require P^ > Aqq]-, and this is the region 

we will consider in details. To this end we split the integral over P^ in Eq. (H) into /p2 + /p2 

' mill 

with some Pq > Aqqj-,. In the first integral we use the VDM-like photon structure function, 
while the whole dependence on P^ is contained in the second one: 
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(10) 







III. QCD CALCULATION OF ELECTRON STRUCTURE FUNCTION 

From the theoretical point of view the QCD behavior of structure functions is most easily 
analyzed in terms of their moments, defined as 

1 

f{n) = Jdxx^-'f{x) (11) 



for any function /. 

The electron structure function we are going to investigate has the form of the second 
term of Eq. ([lOD and its moments read 

ti 

n{n,t,h,to) = f^in) Jdt,rk{n,t,t,), (12) 

to 

where 

q2 p2 p2 p2 

^ = log T2 ' ^0 = log , tl = log -2 , = log (13) 

^^■QCD ^^QCD ^^QCD ^^QCD 

and explicit to argument of /Kn, t, ti, to) is to remind on the dependence on 'auxiliary' scale 

p2 

Our task will be to integrate (n, t, t^) over and to construct master equations for 
the electron structure function. In order to introduce notation and get some understanding 
of the energy scales involved, let us first briefly remind the derivation of the virtual photon 
structure function. 

The master (DGLAP) equations p| read 

dt ~ 2n ^ 

+ ^PMfl-{n,t,t,) + ^P,c(n)/2(n,t,t,), (14a) 

^^^^^^ = ^Po.{n) E [/;(n, t, t,) + f^{n, t, t,) 
ol 

+ —PaG{n)f^{n,t,t,), (14b) 

ZTT 
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where 



Pg^{x) = 3[x'^ + {l-x) 



(15) 



is the photon-quark sphtting function and Pik are the QCD (Altarelh-Parisi) sphtting func- 
tions. In these equations the photon virtuahty (t-y) is fixed and can be thought of as an 
"external" parameter describing the state, QCD content of which depends on t and n. The 
inhomogenous term in Eq.( |14a| ) makes the difference with QCD equations for hadrons. 

The standard method of solving the evolution equations is to decompose first the struc- 
ture functions into singlet and non-singlet components. To simplify the discussion we will 
present formulae for the non-singlet part only. Defining the non-singlet part of a structure 
function as 



1 "* 

fq/q = fq/q ~ TTT' ^^ifq + fq 



(16) 



we obtain 



where 



dt 



27r 



(17) 



1 



(18) 



In the leading log order of QCD 



a{t) _ 2 
~ ]3oi 



with /5o = 11 — 2nf/3 for Uf flavors, and Eq. ([T7|) becomes 

- = dq^in) - -^^^fg/gin, t, t^), 



dt 



where 



and 



dg^{n) 



^{el-{el))Pq,{n) 



(19) 



(20) 



(21) 



The general solution to this differential equation reads 



dqy{n) t 

1 + dqq{n) 
+ fq/q{n,t',t. 



1 - 



t 



l\ l+dqq{n) 



(22) 



(23) 
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As shown for the first time by Witten [|l| , the first term is characteristic feature of the photon 
structure function and for large t it dominates, resulting in linear growth with t. The second 
term in this solution depends on the 'input' measured at some t' and is analogous to the 
leading-log QCD evolution of hadronic structure functions. For low t' and the measured 
structure function fg/q{n,t',t^) can be identified with the VDM-like fgj^{n,t',t^) discussed 
in the previous section. 

Eq.(|23|) gives no explicit prediction on the dependence of f'^^-{n,t,t^) except for the 
fact that the asymptotic {t ^ t') solution is independent of 

Another method of calculating QCD structure functions is the ladder expansion. The 
corresponding diagram for the photon structure function is shown in Fig.|^a. Actual cal- 
culations should be performed in the axial gauge but here we will integrate only over the 
quark emitted by the photon, which is as simple as 

'3-' 2 

fUn,t,t,) = d,,{n) J -^Un,Q',e). (25) 




a) b) 

FIG. 2. General ladder expansion for the photon (a) and the electron structure function (b) 

The QCD 'structure function of a point-like virtual quark' fqq{n,Q'^,k'^) corresponds 
to the ladder diagram of Fig. ^ without the lowest quark rung. The latter is explicitly 
integrated in Eq.(^) with the electromagnetic 7-g coupling contained in dq^{n) and 
coming from the quark propagator. The only difference with pure QCD is that here the 
coupling is electromagnetic and does not depend on k"^. Changing the integration variable 
to r = \og{k^/AQQ^) and using the QCD formula 
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T 



dqq{n) 



(26) 



we arrive at 



i /^\dqq{n) 



dq^{n) t 

1 + dqq{n) 



1 



\l+dqq{n)' 



Formally this result equals to the solution of master equations Eq.(P^ with t' 
0. 

Let us explain the implicit assumptions made in the derivation of Eq.(|2 
the strong ordering of virtualities in the ladder expansion 

P2 < < . . . < Q2 



(27) 
(28) 



and 



Thanks to 



(29) 



we could integrate over the whole range of quark virtualities k"^. Within QCD, however, this 
can be done only if the photon virtuality > Aqqq. Moreover we have used the fact that 
in perturbative calculation such photon has a point-like coupling to quarks. In other words 
the photon of virtuality P^ > Aqq]-, has no QCD structure at the scale = P^. We see, 
thus, that Eq.(l^) is valid for any P^ while Eg . (1281) for P^ > Aq^j-, only. This is exactly the 
reason for introducing the intermediate scale Pq Eq.(|TO|). 

In the following we will assume that P^ is large enough for Eq.(^) to hold. With this 
assumption the integration over P^, as in Eq.(|T^, is straightforward and results in the 
following expression for the electron structure function 



r,{n) J dt,f]/-{n,t,t,) 

to 

f^{n)dgy{n)t 

1 + dqqin) 



t 



2 + dqq{n) 



tl — to 



Any 



(30) 



This result corresponds to the diagram depicted in Fig. 0b, where the range of photon 
virtualities is controlled by imposing a limit on momentum transfer to the outgoing electron. 
Eq.(pOD depends on three scales t > ti > to with to kept fixed. In order to find formulae 
for large t {t ^ to) we have to specify the relation between ti and t. Let us consider two 
extreme cases: ti <^t and ti = t. 



'Exclusive' case: ti <tit 



fq/q{n,t,tl,to) 



qi 



n] 



1 + dqq{n) 



t{ti-to). 



(31) 
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• 'Inclusive' case: ti = t 



fldn, t) ^ r^^.in, t, t, to) - t\ (32) 



where we have dropped the last two arguments of the 'inclusive' structure function. 

To obtain the full result for /^^^ we still have to add the integral over the photon virtual- 
ities below Pq. To this end we use Eg. (p3D with experimental input replaced by a VDM-like 
parametrization at t' = Iq: 



dqy{n) t 



qq 



n 



1 - 



T 



+fiJ^in,t,,t,) (I) . (33) 

As discussed earlier fgJq{Ti,to,tj) should decrease with increasing and vanish for > to 
(see e.g. Q for a phenomenological parametrization). Thus for t ^ to we get the unique 
prediction independent of the 'input' at low to 

f^in) I dt,r,in, t, t,) ^ G^^^^2!^(to - t^i.), (34) 

^min 

where t^m = log(P^.^/A^cD)- 

For large t this low contribution grows linearly with t and can be neglected in the 
'inclusive' case. Eq.(plD becomes now 



This is exactly the Weizsacker- Williams formula Eq.(y) with asymptotic solution Eq.fplD 
used for the photon structure function. 

Much more interesting is the 'inclusive' case Eq.(^). To understand this QCD prediction 
let us look first at the photon structure function. There the effect of QCD evolution can be 
seen by comparing the full result Eq . (pS]) with the Quark Parton Model (QPM) limit. We 
reach this limit by taking Aq^j-, (a(t) — > 0) in Eq. (PBD , which results in 

g2 

fq/qif^^ ^7)Iqpm = dq^{n) (t - t^) = dg^{n) log — (36) 

7 

Thus we see that for large t the net effect of the QCD evolution on the photon structure 
function is to change dq^{n) into dq^{n)/{l + dqq{n)). The dependence on t remains the same 
but the structure functions (transformed back to the a;-space) have different dependence on x. 
As usually the QCD evolution 'shifts' the distribution towards lower x values. Analogously, 
for the electron structure function the QPM limit of Eq. (|30|) reads 
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/^/^-(n, t, ti, to) Iqpm ^ llM^n^^t, - to) (2t - ti - to) 



' ^'^"■"■"log^llog^ + log^J (37) 



The reader can easily check that this resuh corresponds to the integral Jt^dt^ of the QPM 
formula for the photon structure function, Eq.(p^). 

Comparing the QPM result with the QCD formulae Eq.(PB[) and Eq.(|52|) we see that 
the effect of QCD evolution is to multiply the moments of the electron structure function 
by 1/(1 + dqq{n)) in the 'exclusive' case and by 2/(2 + dqq{n)) in the 'inclusive' case. After 
transforming back to the x-space, this means that QPM, 'exclusive' and 'inclusive' electron 
structure functions all have different dependence on x. In particular the 'inclusive' solution 
f^/g{n,t), which corresponds to a standard structure function, gets modified analogously to 
the photon case but by another factor. In this sense the electron acquires an anomalous 
component from QCD. 

So far we have discussed the non-singlet solution. The singlet case goes along the same 
lines but will not be presented here. Instead, we construct in the next section the evolution 
equations which can be solved in the x-space. 



IV. EVOLUTION EQUATIONS 

Let us first show that Eq. (pOD is the general solution to the following master equations 

(38a) 



dfq/q{n,t,tl,to) , , dqq{n) .\ 



dt 

dfq/q{n,t,ti,to) 



dti 



f»Pq/-q{n.tM) 



dqein)t 
1 + dqq{n 



1 



(38b) 



where dqe{n) = f^{n)dq^{n). Note that the second equation is just the derivative of Eq.(|TUp 
with Eq . (|28|) inserted for the photon structure function. A general solution to Eq.( p8a| ) can 
be written as 



rq,q{n,t,h,to) = C{n,t{)t- 



+ 



dqejn) 

1 + dgq{n] 



(ti -to)t. 



(39) 



So far C{n,ti) is an arbitrary function of ti. Nb. if ti remains constant when t — > oo the 
second term of Eq.(p9[) gives the asymptotic solution for the 'exclusive' case. 
From the second equation Eq. (|38b| ) we obtain 



C(n,ti) = C(°)(n) 



dge{n) ti 



dqq{n) + 2 



[1 + dqq{n)] [2 + dqq{n)] 



(40) 



with arbitrary C^^\ 



n] 
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Now the general solution to Eq.(|38D reads 



1 + dgg{n 



2 + dqq{n) 



(41) 



Upon imposing the boundary condition fq/g{n,t,to,to) = we recover the formula Eq.(pOD 
derived in the previous section. 

The complete set of master equations in the x-space analogous to Eq.(^) can be obtained 
by changing the products of moments into convolutions and observing that there is no 
inhomogenous term for the gluonic component — c.f. Eq.(14). Suppressing explicit to 
dependence in the function arguments we have 



where Ap = P„, 



dt 



dn,g{tM) 

dti 
dt 



dti 
t. 



Up 



-Pqe\tl 



tn 



27r 

01 a 



(42a) 
(42b) 



9=1 



a 



+ — PGG®fGit,t,] 

t^fcitM). 



(42c) 
(42d) 



These equations, as they stand, are not very useful because they contain both electron 
and photon structure functions. We will present now how a closed set of equations for the 
electron structure function is formed in the large t region. As already discussed previously, 
the large t limit depends on additional assumptions on ti vs. t dependence. Quite generally 
one can take ti to be some function of t defining a 'path' in the t-ti plane along which the 
large t limit is approached. Denoting ti = g{t) we have 



dt 



A simple choice for g{t) is 



jj. ~ aj. + y {'^) 



dt 

dmt,t,] 

dt 



dmt,h 

dti 



+ g\t)t®fl{tM)- 



g{t) = (1 — a)ii + at 



(43) 



(44) 



where ti is constant and < a < 1. The two extreme cases considered in the previous 
section correspond to following choices for the parameter a 
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'Exclusive' case 

a = 0, i.e. g{t) = ti. As g'{t) = the second term in Eq.(43) vanishes. The solution 
grows hnearly with t {c.f. Eq.(^)). The resulting master equations have a constant 
(t-independent) inhomogenous term 



dt 



dt 



-Pqe [tl 



tn 



2tt 

OL ^ OL 



271 



a 



2n 



Pr 



Gq 



E [fqit 



2n 



9=1 



a 



+^PGG^m,h] 



(45a) 



(45b) 



• 'Inclusive' case. 

a = 1, i.e. g{t) = t and g'{t) = 1. Now the virtual photon structure function f^{t,t) 
vanishes because both arguments are equal. We arrive at the equations which are 
formally the same as in the 'exclusive' case but with ti set to t and tmin neglected. The 
master equations have now the inhomogenous term proportional to t, which results in 
a different x-dependence of the solutions. 

According to Eq.(^) and Eq.(|32D the asymptotic solutions to these equations have the 
form 



27T 



(46) 



with tl constant and ti = t for the 'exclusive' and 'inclusive' case, respectively. Here we 
have taken to = (Pq = Aq^^j^). Substituting the leading-log formula Eq . ([T9|) for a(t) we 
end up with t-independent integral equations 



Pfq/qiz) 

p/g(^) 



Pqeiz) + ^ I- 



X 



Pqq{^) fq/q[ ^ 



X 



dx 



X 



X 



P.Gix)ra[- 



k=q,q,G 

where p = 1 or p = 2 for the 'exclusive' or 'inclusive' case, respectively. 



(47a) 
(47b) 
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z 

FIG. 3. Asymptotic n-quark 'exclusive' and 'inclusive' distributions zf^{z^Q'^)/a^^ for 
Q"^ = lOOGeV^ and = IGeV^ for 'exclusive' case. Calculations have been done for 5 flavors, 
Aqcd = 0.2GeV and = A^^cD- 

We solve these equations numerically using the method described in In Fig.^ we 
show the 'exclusive' and 'inclusive' solutions for the distribution of u quarks in the electron 
at = lOOGeV^. There is still the contribution from the photons of low virtuality (P^ < 
^qcd) which should be added, as it has been done in Eq. (pS]) . At asymptotically large t = 
log((5^/AQ(jj3) it modifies the 'exclusive' solution only, by changing the ti = \og{P^ / Aq^j^) 
factor in Eq.(PB|) into log(P^/Pj^i^). At finite Q^, however, this contribution is non-negligible 
for 'inclusive' case and should also be added. As can be seen from Eg . (|30|) taking the 
intermediate scale Pq = Aqq^) (to = 0) sets the proper normalization for this correction]^ at 
finite t. Thus the both curves in Fig. ^ get shifted by the same amount. 

V. SUMMARY 

In this paper we have analyzed the parton content of the electron within perturbative 
QCD. We have shown that electron acquires an anomalous component from QCD, anal- 
ogously to photon. We have constructed the evolution equations for the 'exclusive' and 
'inclusive' electron structure function. These two cases correspond to 'anti-tagging' and 
'no-tagging' experimental conditions, respectively. The evolution equations can be solved 



^At value considered here the low contribution is of the same order of magnitude as the 
QCD results in Fig.|. 
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numerically in the x-space in the asymptotic region. As an example we have shown the 
u quark distribution inside the electron. 

The results presented here are leading-log QCD solutions valid at asymptotically large 
. At finite the next-to- leading corrections, as well as non-perturbative contributions 
including hadronic component of the real photon will modify the results. Despite these 
inaccuracies it would be interesting to compare these predictions with experiment. On one 
hand, the data for the electron structure function should be much more precise than the 
ones used for the photon structure function. On the other hand, the improvements on 
the theoretical side can be done in a similar way as for the photon structure function — 
higher order perturbative QCD effects, as well as phenomenological parametrizations can be 
plugged into the evolution equations Eq.(|45|). 

The author would like to thank Jerzy Szwed for numerous discussions and for critical 
reading of the manuscript. The hospitality of DESY Theory group, where part of this work 
was done, is also acknowledged. 
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